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Abstract 



Oh! 



The problem of two electrons in a two-dimensional random potential is addressed 
numerically. Specifically, the role of the Coulomb interaction between electrons on 
c/2 \ localization is investigated by writing the Hamiltonian on a localized basis and diago- 

nalizing it exactly. The result of that procedure is discussed in terms of level statistics, 
the expectation value of the electron-electron separation, and a configuration-space in- 
verse participation number. We argue that, in the interacting problem, a localization- 
delocalization crossover in real space does not correspond exactly to a Poisson-Wigner 
crossover in level statistics. 
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PACS 71.23An; 71.30+h; 71.23-k 

§1. Introduction. The problem of two interacting particles (TIP) in a random potential 
has received much attention in the last few years. The focus has been primarily on TIP in one 



dimension (ID), where most investigators have dealt with particles interacting via an on-site 
potential. (For a concise summary of the various approaches used, and results obtained, we 
refer the reader to the articles in ref. 0.) The reason for all the attention (and controversy) 
is the result, first found by Shepelyansky 0, that the TIP actually propagate coherently 
through a length much larger than the one-particle localization radius, which can lead to an 
enhancement of transport ||. We address the related problem of localization in 2D systems, 
and how it is affected by a long-range electron-electron interaction (EEI). 

The more general problem of the interplay between disorder, interaction and quantum 
tunneling, and their combined effect on electronic localization is not new. While the Hartree 
Coulomb repulsion introduces an additional random energy and thus enhances localization, 
the possibility that quantum correlation due to the EEI may act to delocalize the electrons 
was proposed twenty years ago by Pollak and Knotek [Q, and Pollak JJJ, but a firm answer 
has not been achieved yet. Computationally, the main difficulty in the finite-electron-density 
problem is the huge phase space for systems of reasonable size || . Existing work [|, [7], ||, ||, |H| 
resorted to various approximations. In contrast, the TIP problem for reasonably large sys- 
tems can be solved without such approximations: double occupation of sites can be accounted 
for, spin and exchange included, and the entire phase space can be examined. A motivation 
for the problem considered here is the ability to make inferences about approximations made 
in investigations of finite-density and few-electron systems ]TT| , [T2| , [T3"| , |T1| , [T3], |HJ . We hope 



furthermore that the work may contribute to insight into the mechanisms at play in the 
experimental reports on an observed metal-insulator transition (MIT) in 2D [17]. 
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In the present paper, we deal with two electrons in a random potential interacting via 
a long-range Coulomb interaction, and investigate numerically the effect of that EEI on 
electronic localization. A tool used frequently to assess localization has been the distribution 



p(s) of nearest-level spacings s. In the absence of interactions, it has been shown [18| that 
p(s) shifts from Poisson to Wigner as the system goes from being strongly localized to 
delocalized. (The salient difference between the two distributions is that Poisson is maximal 
as s — > 0, while Wigner vanishes there.) For interacting systems, such correspondence has 
never been proven. Still, level spacing statistics has been used often to study localization also 
with interactions. This may be reasonable because interaction is not relevant to the logical 
connections between the localization <-> derealization and Poisson <-» Wigner transitions: 
off-diagonal energies cause hybridization of site functions and thus derealization, while at 
the same time they increase level repulsion and thus eliminate small level spacings. This is 
not to say that the same criteria for the Anderson transition that were established for non- 
interacting systems can be automatically taken over for interacting systems. An interesting 
case in point is a study of the Two-Body Random Interaction Model (TBRIM). Georgeot 
and Shepelyansky [[H| found that in that model a huge number of non-interacting eigenstates 
contribute to the interacting eigenstates (indicating possible derealization in real space) even 
in situations when p(s) is close to Poisson. On the other hand, Jacquod and Shepelyansky 
20| used, for the same model, the Wigner-Poisson transition to establish a transition from 
integrability to chaos. 

A quantity used commonly to study derealization in non-interacting systems is the 



inverse participation number P l ||21|| , which measures over how many sites the (one-particle) 



wavefunction spreads. Here, where we allow for interactions, we use an analogous quantity, 
the configuration-space inverse participation number R (see ref. |22] and eq. @ below) 



which measures over how many configurations the many-particle wavefunction spreads. In 
the limits of strong localization and of complete derealization it is easy to see the connection 
between P~ x and R: the strong-localization limit, R = 1, implies the presence of a single 
configuration, which in turn implies that each particle is localized on a single site. In the 
complete derealization limit, R = 1/N (where N is the total number of configurations), 
the wavefunction extends uniformly over all configurations, implying that the particles are 
uniformly spread over the system in real space. 

Previous studies on 2D random systems with interactions most relevant to this work 



include finite-size scaling of three and four [l 1| , and two |TJL O, O, [15|, [16] spinless elec 



trons. All these works have concluded that the interaction enhances derealization. In [III 



crossover from Poisson to Wigner was found, while [jl^, [L3|, [0| reported a sharp transition 
with an identifiable critical point. We shall see that the enhanced accuracy of this study 
does not alter the conclusion that interaction enhances derealization through most of the 
energy domain. 

In the next two sections, we explain the details of our approach, and describe the differ- 
ence with other authors' methods. In the fourth section, we present our results, and then 
we offer in the conclusion some final remarks. 



§2. The Model. In the standard approach (see, e.g. refs. |TT| , |T5|| ) to the problem 
at hand (sometimes referred to as the Quantum Coulomb Glass 0), one considers spinless 
electrons on sites of a lattice and represents the Hamiltonian in the basis of local wavefunc- 
tions: 

Hqcg = t^2(at a j + h - c -) + X^ a l a * + U^2(a\aia]aj + a\a,ja)a,i)/\ri - fj\ (1) 

{ij} i i>j 

where {ij} denotes nearest-neighbor (n-n) sites, the operator a\ (a*) creates (destroys) a 
(spinless) electron at site i; 6j is a random site energy chosen in the range — W/2 < 6{ < W/2 
with uniform probability; U is the strength of the Coulomb interaction; r is a position vector. 
Both U and W are taken as independent parameters. 

Our model differs primarily in that we include spin, and treat more extensively the 
overlaps. The latter is done in two ways: we do not neglect fluctuations in the n-n overlap 
integrals (due to the differing charge environments at different pairs of sites), and include 
overlaps over other than just nearest-neighbors (as we cross the MIT from the insulating 
side, the increase in the number of important distant sites may be more decisive than the 
decrease in the n-n overlap.) 

To accomplish the above we write the two-electron wavefunctions in terms of a basis set 
of appropriately symmetrized products of the one-particle local wavefunctions 

$ y = 0jf xa (A) ,and (2) 

*« = 4 A) x^ } > (3) 

where i, j are site labels, 0y (a) is the spatial (spin) part of the two-electron configuration, 
and the superscripts S and A indicate the symmetry of the wavefunction under particle 
exchange. The <f>ij are constructed in the usual manner, i.e. by symmetrizing (or anti- 
symmetrizing) products such as ipi(l)ipj(2) of one-electron orbitals (centered on i and j) ip 

(S) (S) 

for electrons 1 and 2. Double occupation of sites comes in through (fr^ . Clearly, the 4>\j 

[cplj ) are singlet (triplet) configurations. We take (p ~ exp(— r/as), with r the electron-core 
distance and as the microscopic (Bohr) radius corresponding to those orbitals. 
We write the Hamiltonian 

2 2 

H = J2(T a + V a + e a ) + —,with (4) 

T a = — Vl, 

2m* 

Us 1 

V a = -e-q s J2 > and 
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where a labels the electrons, T, V, and e are the operators corresponding to the kinetic 
energy, the interaction of the electrons with the cores, and the random potential, respectively; 
e is the electronic charge; k is the dielectric constant; t\i is the electron-electron (e-e) 
distance; m* is the electron mass; q s is the charge on a site; n s is the number of sites in 
the system; e« is a random energy chosen as in eq. ([!]), with W equal to the n-n Coulomb 
energy ec (i.e. we set W = U = 1 in units of a n-n Coulomb energy). We vary the site 
concentration by changing the n-n distance r nn as a parameter. 

Analytic solutions for (<%- \H\(jr kl ) and {4>\j \H\<p kl ) were derived and written in terms 
of one-, two-, three-, and four-center integrals. The procedure is very tedious but straight- 
forward. The equations are not given here for reason of space, and because the TIP case 
is not general; the equivalent expressions for systems with an arbitrary number of spinless 
fermions were published in [f2I|. The integrals were performed numerically for processes 
{ij} — > {kl} which involved (i) no electron transfers (i.e. diagonal matrix elements); or (ii) 
a one-electron n-n transfer; or (Hi) a next n-n transfer; or (iv) a next-to-next n-n transfer; 
or (v) two simultaneous n-n transfers. All other off-diagonal matrix elements were set to 
zero. H does not depend on spin, so a enters into the picture only by dictating the symme- 
try (under particle exchange) of the corresponding 0y. Naturally, the off-diagonal matrix 
elements between faj of different symmetry are set to zero; thus the matrix H splits into two 
blocks which we diagonalize separately using standard techniques. The result is the set of 
eigenenergies Ej, and the corresponding states $/ in which the 0y come in with amplitudes 

if' = £ 4M 5) . »a *<"» = £ 4S4 A) . (5) 

{ij} {M 

For simplicity we refer below to &} as "singlets" and $]- as "triplets". 

In this work we take W to always equal ec; thus, W = ec ~ l/r nn . Since r nn is a 
parameter here, W changes accordingly. We do this because then the effect of interaction is 
most important, and also because this condition prevails in many experiments, for example, 
in the vast experimental literature on impurity conduction at moderate compensation [^3|. 
In impurity conduction W corresponds to the Coulomb energy over the mean n-n majority- 
minority ion distance, while U corresponds to the Coulomb energy over the mean carrier- 
carrier separation. The two are quite similar unless the compensation is either very small or 
very large. 

The present model also differs from other models in how it deals with electrical neutrality. 
In most works neutrality is achieved by placing a compensating charge q s = +\e\n e /n s on 
every site. However, since the microscopic radius of the electrons is usually determined by 
the charge q s of the core, this radius would strongly depend on how electrical neutrality is 
established. (For the TIP problem n e is fixed, but n s is generally taken as a parameter; 
thus, different models for compensation may in fact yield different results.) To avoid such 
problems, this work does not consider the compensating charge explicitly. To make contact 
with experimental situations, we place a core charge of magnitude +|e| at each site. The 
macroscopic radius is then determined by an effective mass and a dielectric constant. Since 



there are more sites than carriers, the system is not neutral. A question we need to answer 
is whether the properties of interest here, i.e. the degree of localization and its dependence 
on various parameters, are fundamentally affected by the lack of neutrality To answer 
this question we consider whether one can construct some model for the compensating 
charge such that it has no effect on the results of the computations as far as the localization 
properties are concerned. The answer is that such a model exists, namely a spatially uniform 
distribution of the compensating charge (which may be placed on a separate plane, not unlike 
the situation in gated devices). This charge merely contributes two constant energies: a self- 
interaction, and the interaction with the two electrons. We may thus conclude that the lack 
of neutrality per se does not affect our results. However, it should be borne in mind that 
in reality the particular way in which charge neutrality is realized may affect somewhat the 
localization properties. 

§3. Procedure. We set up "samples" on the computer with n s (= L x L) sites arranged 
on a lattice, and diagonalize eq. (|4]) for the parameters L and r nn . For definiteness, we 
take as and k to be 10A and 3 respectively, which yields and effective mass m* = 0.16m, 
with m the electron mass - these values seem appropriate for 2D systems. Cyclic boundary 
conditions are used. In what follows, we use a# as the unit of distance. In every case the 
number of samples was sufficient to obtain at least 1.5 x 10 4 levels for each pair {n s , r nn }. 

As one measure of localization, we compute p(s). (We dropped up to 100 states from 
the band edges.) As in flT^j , a parameter r\ = (var(s) — 0.273)/(1.0 — 0.273) is computed as 
a measure of how close p(s) is to a Poisson (77 = 1) or Wigner (rj = 0) distribution. As an 
alternate measure of localization we compute the configuration-space inverse participation 



number as in ref. 22 



^ = £1^/. (6) 

In addition to p(s) and Ri we examine the behavior of the quantity 

^i = ^2 r ij\A I , ij \ 2 , (7) 

m 

which is the expectation value of the e-e separation when the system is in state $/. This 
expectation value is computed here for a direct glimpse at the behavior of the $/ in real 
space. One expects that in the localized regime A/ is strongly correlated with Ej due to the 
EEI (larger A/ correspond to smaller Ej), whereas A/ should become essentially independent 
of Ej as configuration-mixing increases. In the absence of interactions, such correlations 
should of course not be present. 

To measure the importance of secondary tunneling processes, we investigate the following 
quantities: (i) the width w, and the average t nn of the distribution of off-diagonal elements 
corresponding to n-n processes; (ii) the sum S^ of all off-diagonal elements corresponding 
to n-n processes; (Hi) a similar sum S^ 2 ) for all next n-n processes; and (iv) a sum S^ 
of such matrix elements which correspond to either next-to-next n-n or two simultaneous 



one-electron n-n transfers. We note that the sums S^'"- 1 account in a crude way for not 
only the magnitude of the off-diagonal matrix elements, but also for the size of the phase 
space occupied by those processes {e.g. even though a matrix element corresponding to 
a n-n transition may be significantly larger than one corresponding to two simultaneous 
n-n transfers, the number of matrix elements of the latter type is much larger - and thus 
such processes may contribute significantly to coherent tunneling). Also, we note that the 
standard approach assumes t nn 3> w, S^ ^> S^^ and S^ ^> S^ 2 \ 

§4. Results and discussion. Double occupation of sites does not, of course, occur 
in the triplets case, and a single band of eigenenergies results from the diagonalization 
procedure. The width of this band naturally increases with decreasing r nn because the 
configurations hybridize, and level repulsion increases. For singlets, however, at large r nn two 
such bands separated by a gap result: a lower band arises from hybridization of configurations 
in which sites are singly occupied (to which we will refer as type-s configurations) - as in the 
triplets case, and an upper band arises from configurations in which sites are doubly-occupied 
(type-d configurations), which for the most part do not hybridize; with decreasing r nn , the 
gap narrows as (mostly type-s) configurations hybridize, and level repulsion increases. As 
r nn is decreased further, type-d configurations start to mix with each other, and with type-s 
configurations, until eventually (for r nn <~ 4), the gap actually disappears, and the two 
bands start to join together. We present results here for r nn = 4-12. In what follows, except 
for our discussion of Xj, our results refer to the triplets band, and to the singlets' lower band. 

We first present the basic conclusions regarding the effect of the EEI on localization. 
Fig. [I] shows a comparison of i](r nn ) between the interacting and non- interacting systems 
[i. e. eq. (|j) without the last term] . It is very clear that this criterion shows EEI to enhance 
derealization - while with EEI derealization occurs at r nn ~ 10, without EEI derealization 
is out of the range of the figure, namely somewhere at r nn < 7. We note that there is no 
clear small-size scaling behavior, suggesting a crossover (r nn ~ 9-11), rather than a critical 
transition. This is in agreement with [11] , but differs from [[[2, [13], [15], 16|. The differences 



may be due to differences in models and choice of parameters 04] . 

We present in figs. |] and our results for Rj with L = 7, but we checked that the plots for 
other values of L are qualitatively the same. To evaluate the importance of spin, we plot in 
fig. [| the number Ri for singlets and triplets. In both cases, the energies are measured from 
the bottom of the band, and have been normalized by the energy difference between the top 
and the bottom of the band. The most noticeable features in fig. are: 1) in the localized 
regime [fig. H(a)] there is no discernible difference between singlets and triplets; 2) For r nn 
close to the transition and into the delocalized regime, the lowest energy states of the singlets 
are always more delocalized than the lowest energy triplets. This is exemplified by figs. @(b) 
and 0(c). This difference between the singlets and the triplets seems to be larger the smaller 
the r nn . 3) In the well-delocalized regime there appear intriguing fluctuations of Ri with 
energy, for both singlets and triplets, as exemplified by fig. 0(c). The fluctuations are not 
random - there is a definite correlation in Rj over a finite distance in energy. This persists 



from one realization of random energies to another, and persists also for increased L, though 
with a different characteristic correlation width. Possibly the fluctuations are connected with 
interference effects that come into play, due to our use of cyclic boundary conditions, once 
the extent of the wavefunction spans the entire system. Interestingly, similar fluctuations 
were observed p5j in a real-space inverse participation ratio for TIP in the Harper model 



where the non-interacting eigenstates are delocalized. 

In fig H p is the average (over singlets and triplets) Rj with EEI to the average Rj 
without the interactions, i.e. Ri is computed with and without the last term in eq. (^), 
and the results are averaged over singlets and triplets - p is given by the ratio of the two 
averages. Thus, p < 1 indicates derealization (in configuration space) induced by the EEI, 
whereas p > 1 points to interaction-induced localization. Over most of the regime of Ej the 
plots confirm the previous result about the effect of EEI on derealization. However, at the 
extremes of the energy regime, the opposite appears to happen - the wavefunctions for the 
interacting system are more localized, at least in configuration space. (It has been shown 
26| that EEI can, under certain conditions enhance either localization or derealization, 



and this in fact has been found to be the case [lQfl.) We discuss the enhanced localization 



at the band edges in the context of f2q| . The enhancement or supression of localization 
depends on the general effect of the EEI on the ratios Hjj/Du, where Hjj = (&i\H |$j) 
and Djj = (&i\H \§i) — (&j\H |$j). [Configuration indices J, J here replace the pair indices 
i, j of eqns. (2) and (3).] If the EEI enhances Hu/Du it supresses localization, in the 
opposite case it enhances it. Hjj and Djj are both separately enhanced by the EEI 



so 



the effect on localization depends on which is enhanced more. Now consider, for example, 
the lower band edge, i.e. the very low-energy end of the spectrum. The level spacing, and 
thus Dij, is enhanced by the interaction because it excludes configurations with nearby 
electrons. The enhancement of Hjj by interactions comes to a large degree from correlated 
n-particle (here 2-particle) transitions, but this effect is important mainly when the electrons 
are reasonably close-by A somewhat parallel argument can be made for the upper band edge. 
It is interesting to note that similar effects of enhanced localization were reported in [ES[ for 
rather different systems. 

We consider next the question of e-e separation due to interaction as measured by A/. 
Without the EEI the positions of the two electrons are uncorrelated so A/ is expected to vary 
at random from configuration to configuration and we ascertained that this is indeed the 
case. Fig. |] shows polynomial fits for Xj(Ej) with EEI for L = 7. (Again, we checked that 
the results for other values of L are qualitatively the same.) Fig. [|(a) corresponds to singlets, 
and fig. |](b) to triplets. The numbers on the graph refer to the values of r nn . In the singlets 
case we now include both bands. This figure is similar to results presented in ||27|| , but here 
we present an average over many realizations of the random potential. In figs. |](a) and (b) 
the energies are measured from the bottom of the lower band. In fig. |](a) the energies are 
normalized to the energy difference between the top of the upper band the the bottom of 
the lower band; for comparisons between figs. |](a) and (b), the triplets' energies have been 
normalized so that the top of their band coincides with the top of the singlets' lower band. 



The results shown in the figure can be conveniently discussed by dividing the plots into three 
groups: 1) The singlets [fig. f|(a)] include states at the high energy end of the figure which 
do not have a triplet counterpart [fig. (|(b)]. As discussed above, these states correspond to 
double occupation of sites: for large r nn , A/ is quite small and the energy large, as expected 
for two electrons residing on the same site. As r nn decreases, some states still have the large 
energy, but now Xj is larger (as is evident from the steep lines with r nn = 4,5). These are 
likely to be states with a reasonably large component of a doubly-occupied site, hybridized 
with configurations where electrons are more remote from each other. The other features of 
fig. |] are common to the triplets and the singlets and correspond to states where double 
occupation of sites is minimal or non-existent. 2) The curves representing large values of r nn 
(small overlaps) are quite steep, i.e. A/ decreases sharply with increasing Ej. This is easily 
understandable as an increase in the repulsion energy with a decreasing distance between 
occupied sites, i.e. decreasing A/. 3) The dependence of A/ on Ei weakens as r nn decreases 
and configuration-mixing increases - A/ becomes nearly independent of Ej and quite large 
for r nn = 5. The crossover takes place around r nn = 6-9 in both fig. f|(a) and (b). 

We note that the crossover in \j is lower than in fig. |l| In this model, W/t = e Tnn /[r nn (l + 
r nn )}, so a small difference in r nn implies a comparatively larger one in W/t. This gives 
W/t ~ 3, 5, 10, 20, 40, 90, 200, 500, 1000 for r nn = 4, 5, 6, 7, 8, 9, 10, 11, 12, respectively. One 
must be careful, however, when comparing our results with those of other authors because 
here W/t = e c /t. 

In fig. |5| we present our analysis of the importance of fluctuations in n-n hopping en- 
ergies, and of secondary hopping processes. First we note that the relative width of the 
distribution of n-n hopping energies remains roughly constant with r nn , and decreases with 
system size. As the size of the phase space increases with n s , the relative importance of these 
fluctuations decreases, and probably disappears in the thermodynamic limit. Furthermore, 
secondary hopping processes are not very important in the localized regime; however, they 
are important (as expected) in the delocalized regime. These processes become significant for 
r nn <~ 6, and in fact, next-to-next n-n and simultaneous n-n processes are more important 
than next n-n in the delocalized regime. We surmise that these secondary processes "push" 
the transition towards higher values of r nn . 

We now turn to the specific aspects learned from those features of our model which go 
beyond other works [^, [TB|, [14], |T5[ , namely the inclusion of spin and the more general con- 
sideration of elastic tunneling. We note that: 1) As expected, spin plays no role deep in the 
localized regime (since the exchange energy is proportional to ((pi\(fj) 2 ). More unexpectedly, 
spin also seems to have little importance near the crossover to derealization. Perhaps spin 
plays a larger role when one considers more than two electrons. 2) Type-d configurations 
play little or no role near the crossover - they become important only for r nn <~ 4. 3) 
In the localized extreme (i.e. r nn > 12), the rather short span in energy of the eigenstates 
comes from the random and Coulomb energies; as r nn decreases, the broadening of the range 
of eigenenergies is attributable to the growing off-diagonal energies of H. 4) The standard 
model underestimates the delocalizing effect of the EEI for small r nn especially, as shown 



in fig. [|. Many-electron coherent processes are probably important in the finite-electron- 
density problem. This phenomenon was first pointed out in ||. Secondary processes are 
important especially in the delocalized regime - one should consider them if investigating 
systems around the MIT. 

§5. Conclusions. The model used here yields zero density in the thermodynamic limit, 
and so no definite claims or comparisons with experiments can be made. However, if we 
take the electronic localization-delocalization transition at r nn = 6, this simple model yields 
the critical concentration n c = 10 11 cm -2 . Whereas our systems are quite different from the 



experimental ones (see, e.g. |[28| , where a 2D MIT was reported for high- mobility systems 
at low temperatures), it is interesting to note that the experimental value of n c is the same 
as here. This of course might be just fortuitous. 

Where collective hopping of the two electrons is coherent (r s <~ 7), A can be interpreted 
as a coherence length. It is of interest to note that a crossover of p(s) from Poisson to 
Wigner occurs at a somewhat larger r nn than the crossover from a large variation in Xi(Ei) 
to Xi(Ei) ~ const. We have included in our level spacing analysis most of the eigenstates of 
eq. (f|), as opposed to only states in the middle of the band, which is the customary procedure 
(see e.g. [[II], [12], 0). Since we are including states which tend to be more localized than 
those in the middle of the band (as is evident from fig. §), the effect is that our level spacing 
statistics picture is skewed towards the localized regime. The disagreement between the 
crossover on p(s) and that in A/ would be stronger if we were to follow standard practice. 
We believe that (unlike in the non-interacting case) for interacting systems, derealization 
in real space requires a somewhat larger overlap for n-n sites than does the transition to 



Wigner statistics. This is in keeping with previous work [E2], Z9|, where it was observed 
that the wavefunctions are "swiss cheese-like" (i.e. not compact in real space) without 
EEI, but space-filling with the interactions p0| ; therefore, while the EEI may make the 
wavefunction extend over more sites, it does not similarly increase its spatial extent - A/ 
might require a larger overlap than p(s) (for a crossover to take place) because the EEI 
makes the wavefunction more compact in real space. In a sense, there is a similarity here 



with |y|, where it was shown that interacting eigenstates may contain contributions from 
many non-interacting eigenstates, even when p(s) is close to Poisson. The difference is that 
we get p(s) close to Wigner in some situations when one might expect Poisson. For example, 
from fig. |], at r nn = 9 one would expect p(s) similar to Poisson since Xj depends strongly on 
Ej] however, fig. [I] reveals that p(s) is nearly Wigner. We attribute the discrepancy between 
this work and the result in [|I1| to the different models used. 

It is clear that p(s) is Poisson well inside the localized regime, and Wigner well into 
the delocalized regime; however, the crossover in level statistics and its relationship to a 
localization-delocalization transition in real space is not as well established as it is in the 
non-interacting problem, and requires further study. Furthermore, whereas it is clear that 
deep localization in configuration space implies deep localization in real space, and the same 
holds true for well-delocalized systems, the actual transition need not occur simultaneously 
in the two domains. 
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Figure 1: Distribution p(s) of nearest level spacing. The lines without symbols correspond to 
the non-interacting Hamiltonian. The lines with symbols correspond to the full Hamiltonian 
[eq.g. 
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Figure 2: Inverse participation number in configuration space R with L = 7 for singlets (o) 
and triplets (x), and (a) r nn = 12; (b) r nn = 8; (c) r nn = 6. The energy scale is described 
in the text. 
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Figure 3: The ratio p of the configuration-space inverse participation number R with EEI 
to R without. The numbers on the graph indicate the corresponding values of r nn . p = 1 
for r nn = 12 (not shown). The energy scale is described in the text. 
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Figure 4: Fits for X(E) for (a) singlets, and (b) triplets with L = 7. The numbers on the 
graph correspond to the values of r nn . In (a) the upper bands coincide for r nn = 7-12; in (b) 
the plots for r nn = 11, 12 coincide. The energy scale is described in the text. 
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Figure 5: An evaluation of the importance of secondary processes relative to n-n tunneling. 
The solid lines are plots of w/t nn . The dashed line is S^^/S^ vs. r nn (the same plot for 

(also the same plot for L = 5,7, 10). 



VS. TV 



L = 5,7, 10), and the dotted line is S (2) /S (1 
The vertical axis is given as a percent in all cases. 
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